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The Polyakov-quark-meson (PQM) model, which combines chiral as well as deconfinement aspects 
of strongly interacting matter is introduced for three light quark flavors. An analysis of the chiral and 
deconfinement phase transition of the model and its thermodynamics at finite temperatures is given. 
Three different forms of the effective Polyakov loop potential are considered. The findings of the 
2 + 1 flavor model investigations are confronted to corresponding recent QCD lattice simulations of 
the RBC-Bielefeld, HotQCD and Wuppertal-Budapest collaborations. The influence of the heavier 
quark masses, which are used in the lattice calculations, is taken into account. In the transition 
region the bulk thermodynamics of the PQM model agrees well with the lattice data. 
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I. INTRODUCTION 

While the underlying microscopic theory of the 
strong interaction. Quantum Chromodynamics (QCD), 
is known for several decades, our knowledge on its phase 
structure under extreme conditions is still rather lim- 
ited [l], Q- Understanding its critical properties and 
phase transitions is a major focus which is addressed 
in both theoretical and experimental studies. Heavy- 
ion collisions in such experiments at the highest cur- 
rently available energies are surprisingly well-described 
by ideal relativistic hydrodynamics (see e.g. 0] and ref- 
erences therein). This demands a detailed knowledge of 
the bulk thermodynamics and the equation of state of 
strongly-interacting matter. 

Different regimes of the QCD phase diagram can be 
examined by employing various theoretical methods. For 
example, lattice QCD simulations (see e.g. are appli- 
cable to the finite temperature regime of the phase dia- 
gram but at finite chemical potentials the fermion sign 
problem is still a considerable obstacle. At finite tem- 
peratures recent lattice calculations describe the QCD 
thermodynamics reliably since larger volumes are used 
and the quark masses are closer to their physical val- 
ues. Although most lattice simulations for (2 -I- l)-flavor 
QCD agree that chiral symmetry is restored by a smooth 
crossover transition (Bj there is still an ongoing discussion 
concerning the (pseudo)critical temperatures of the chi- 
ral and deconfinement transition Q. One group (RBC- 
Bielefeld, HotQCD) sees a coincidence of both transi- 
tions 0, Q while the Wuppertal-Budapest group finds a 
larger deconfinement temperature compared to the chiral 
one by examining the peaks of the chiral and Polyakov- 
loop susceptibilities 0, • 

Effective NJL-type models can address chiral aspects 
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of the QCD phase diagram and do not suffer from limita- 
tions at finite chemical potentials. It is expected that the 
critical behavior of such models is governed by the same 
universality class as QCD [ll|. But, due to the lack of 
confinement, they fail in describing the thermodynamics 
of strongly interacting matter, in particular the decon- 
finement transition. However, both transitions become 
accessible with extended NJL-type models which are aug- 
mented with the Polyakov loop. Recently, most studies 
in these Polyakov-loop NJL-type models (PNJL) focus 
on two quark flavors with an effective polynomial or log- 
arithmic Polyakov loop potential |12h21| | . The thermody- 
namics of these models has often been compared to older 
two-flavor lattice data with large quark masses (e.g. [22| ) 
while in the model calculations physical values have been 
used. Recent (2 4- l)-flavor studies in the PNJL [H-i^l 
or Polyakov-quark-meson (PQM) models [28l - [30j indeed 
agree better with the corresponding lattice data although 
the used quark masses still differ. 

Much progress has been achieved in lattice studies for 
(2+l)-flavor QCD thermodynamics, e.g. 0, i, [IS . 
In this work a (2 + l)-flavor PQM model study is com- 
pared with the most recent lattice data of the HotQCD 
and RBC-Bielefeld Isl. Isp and the Wuppertal-Budapest 
(WB) collaborations [10[. In the lattice simulations of 
the HotQCD and RBC-Bielefeld collaborations a ratio of 
the physical strange quark mass to the light one of 10 is 
typically used which yields too large light quark masses. 
These quark masses entail a pion mass of TOtt ^ 220 McV 
and a kaon mass of niK ~ 500 MeV using improved stag- 
gered fermions (p4 and asqtad actions). They observe a 
coincidence of the chiral and the deconfinement transi- 
tion around temperatures T ~ 185 — 195 MeV. Coincid- 
ing temperatures result also in a two fiavor QCD analysis 
if quantum fluctuations arc incorporated, see e.g. [ssHSGf . 

For the chiral order parameter and the strange quark 
number susceptibility we also compare our results with 
lattice data from the Wuppertal-Budapest collaboration 
which they have obtained for physical quark masses on 
Nr = 10 lattices [l^l- In contrast to the other lattice col- 
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laborations they find a chiral transition at temperatures 
lower than the one for the deconfinement transition, i.e., 
Tx ~ 146 MeV < - 170 McV. 

The outhne of this work is as follows: after intro- 
ducing the PQM model for Nf = 2 + 1 quark fla- 
vors, the implementation of various effective Polyakov- 
loop potentials is discussed. The thermodynamic po- 
tential of the PQM model in mean-field approximation 
is derived in Sec. IIIII The resulting chiral and confine- 
ment/deconfinement phase structure and its dependence 
on the model parameters and various Polyakov-loop po- 
tentials are presented in Sec. IIVI A detailed comparison 
for several thermodynamic quantities such as the equa- 
tion of state and various susceptibilities with lattice data 
is carried out in Sec. |Vl Here, the effect of the heavier 
pion masses and different Polyakov loop potential choices 
is explicitly considered. Finally, conclusions are drawn in 
Sec. ED 



II. POLYAKOV-QUARK-MESON MODEL 

In order to address the linkage between the chi- 
ral symmetry restoration and aspects of the con- 
finement /deconfinement transition we introduce the 
Polyakov-quark-meson (PQM) model for three quark fla- 
vors. It is a straightforward generalization of the two 
quark flavor model |36l and is a combination of the chi- 
ral linear cr-model [37| with the Polyakov loop ^(x), the 
thermal expectation value of a color traced Wilson loop 
in the temporal direction 
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{iicV(x)) 



(1) 



The matrix- valued Polyakov loop operator Vix) in the 
fundamental representation of the SU [Nc) gauge group 
contains the temporal vector field Aq, 



V{x) = Pexp 




(ir Ao(a;, r) 



(2) 



where P denotes path ordering and /3 ~ 1/T the inverse 
temperature [ssl ] . Since the Polyakov loop in this model 
is used as a classical variable, implementation details in 
Eq. ([2]) are not important in the present work. 

The expectation value Eq. ^ serves as an order pa- 
rameter of the center symmetry in the heavy quark mass 
limit: it is finite when the center symmetry is sponta- 
neously broken and vanishes for center-symmetric states. 
Hence, $ is finite at high temperatures corresponding 
to the deconfined plasma phase and it vanishes at low 
temperatures in the confined, center-symmetric phase. 
However, in a system with dynamical quarks the cen- 
ter symmetry is always broken but $ still seems to be 
a good indicator of the confinement/deconfinement tran- 
sition. Therefore, we regard $ as an approximate order 
parameter for this phase transition. 



The PQM Lagrangian consists of a quark-meson con- 
tribution and a Polyakov-loop potential Z^(<I>[A], 
which depends on the Polyakov-loop variable 3> and its 
hermitian conjugate $ = (ticV^) /Nc- The coupling 
of the uniform temporal background gauge field to the 
quarks is achieved by replacing the standard derivative 

in the quark-meson contribution by a covariant deriva- 
tive 



D., 



da - iAu 



(3) 



where the SU{Nc) gauge coupling gs is included in = 
gsA'^X°'/2. The usual Gell-Mann matrices are denoted by 
A° with a = 1, • . • , — 1- This leads to the Lagrangian 

CpQM = q{ip-9^5)q + C,n-UmA],^A]) , (4) 

where the column vector q = (u, d, s) denotes the quark 
field for Nf = 3 flavors and A'c = 3 color degrees of 
freedom. For three quark flavors the interaction between 
the quarks and the meson nonets is implemented by a 
flavor-blind Yukawa coupling g and the meson matrix 



05 = Ta (cTa + hsT^a) 



(5) 



where the Ta = Aa/2, a = 0, . . . , 8 are the nine generators 

of the U{3) symmetry with Aq = 1. The generators 

Ta are normalized to TilTaTi,) ~ and obey the 

C/(3) algebra [Ta,Tb] = ifabcTc and {Ta,Tb} = dabcTc 
respectively with the corresponding standard symmetric 
dabc and antisymmetric fabc structure constants of the 
SU{3) group and 



(6) 



fabO — , dabO — \ j -^^ab 



The nine scalar {J^ = 0+) mesons are labeled by the <Ja 
flelds and accordingly the nine pseudoscalar (J^ = 0~) 
mesons by the tt^ flelds. 

The remaining, purely mesonic contribution reads 



A„ - Tr - m2 Tr((/)t0) _ Ai [Tr( 

-A2 Tr (0^)2 + c (det(0) + dct(0^)) 

+ Tr[i7(0 + 0t)] ^ (7) 

with the abbreviation 



= Ta4>a = Ta {Oa + ITTa) , 



(8) 



representing a complex (3 x 3)-matrix for Nf = 3. 

Chiral symmetry is broken explicitly by the last term 
in Eq. ([7]) where 



H = Taha 



(9) 



is a again (3x3)-matrix with nine external parameters ha- 
The [/(l)A-symmetry is explicitly broken by the 't Hooft 
determinant term with a constant strength c. Further 
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details concerning the three-flav or q uark-meson part of with the temperature-dependent coeflicients 
the PQM model can be found in [39| . 

The final form of the effective gluon field potential U 
in Eq. Q is constructed in terms of the Polyakov loop 
variables $ and $. It preserves the center symmetry 
of the pure Yang-Mihs (YM) theory [13, El. However, and 
several explicit choices of the Polyakov loop potential are 
possible which we collect in the following. 



a{T) ao + ai ( y) + f Y 



b{T) - 63 



To 



T 



(13) 



(14) 



A. Polyakov loop potentials 

The functional form of the Polyakov loop potential is 
motivated by the underlying QCD symmetries in the pure 
gauge limit, i.e. for infinitely heavy quarks. The differ- 
ent parameterizations reproduce a first-order transition 
at T ~ 270 MeV in the pure gauge sector for Nc ~ 2> 

The simplest choice is based on a Ginzburg-Landau 
ansatz [Tsl. l36j|. The underlying Z{?)) center symmetry is 
broken spontaneously in the pure Yang-Mills case which 
should be respected by the ansatz. Hence, an expansion 
in terms of the order parameter results in 



^poly 



with a temperature-dependent coefficient 
HT) = a, + aA^)+aJ^)\aJ^)' . (11) 



The cubic $ terms in Eq. ([TU| are required to break the 
C/(l) symmetry of the remaining terms down to Zii). 
The potential parameters are adjusted to the pure gauge 
lattice data such that the equation of state and the 
Polyakov loop expectation values are reproduced. This 
yields the following set of parameters [ij] : 

Co = 6.75, ai = -1.95, 03 = 2.625 , 03 = -7.44 
63 = 0.75 , &4 = 7.5 . 

The parameter Tq = 270 MeV corresponds to the transi- 
tion temperature in the pure YM theory. The extension 
of the potential Eq. to finite chemical potential and 
the remaining ambiguity is discussed in [36|. At van- 
ishing chemical potential the constraint $ = $1' is still 
valid and the potential corresponds to a polynomial ex- 
pansion in $ and $ , including the product term <i><i> as 
in Refs. [MSI. 

An improved ansatz [TJl results in 



rp4 



-ia(T)$$ 



b{T) In [1 - 6$$ + 4 + - 3 ($$)^] (12) 



Here, the logarithmic form constrains $ and $ to values 
smaller than one, in contrast to the polynomial potential 
ansatz. 

Again, the parameters of Eq. (1121) are fitted to the pure 
gauge lattice data, resulting in [ij] 



ao = 3.51 , ai = -2.47 



02 



15.2 



'3 = -1.75 
(15) 

Both fits reproduce equally well the equation of state and 
the Polyakov loop expectation value. 

A third version, proposed by Fukushima [25| . is in- 
spired by a strong-coupling analysis 



'uku 



rp4 



b_ 

Infl 



54e""/^$$ 
-6$$-3($l>)^ 



,4($3 + $3))] (16) 



and has only two parameters a and b. The first term is 
reminiscent of the nearest-neighbor interaction and the 
logarithmic term is again due to the Haar measure as 
in (jl2p . The parameter a determines the deconfinement 
transition, i.e., the transition temperature in pure gauge 
theory, while b controls the mixing of the chiral and the 
deconfinement transition. In this version, the parameters 
are not fitted to lattice data but they reproduce a first- 
order transition at Tq ~ 270 MeV in the pure gauge sector 
for 



664 MeV, &= (196.2 Me V)^ 



(17) 



In a mean-field PNJL analysis with this potential ansatz 
a coincidence of the chiral and deconfinement transitions 
at Tc ^ 200 MeV for vanishing chemical potential is 
found 



25| 



Another significant difference of the last potential 
ansatz is visible in the thermodynamic pressure at high 
temperatures. In this case, both the unconfined trans- 
verse gluons as well as the Polyakov loop which cor- 
responds to longitudinal gluons contribute to the pres- 
sure. But since the Polyakov loop potentials are fit- 
ted to pure gauge lattice data, they thus reproduce the 
pressure of both the longitudinal and the transverse glu- 
ons and might therefore overcount the degrees of free- 
dom in the chirally symmetric phase. However, the po- 
tential ansatz by Fukushima excludes these transverse 
gluon contributions at high temperatures. This differ- 
ence will become important when thermodynamic quan- 
tities, obtained with various Polyakov loop potentials, 
are compared to lattice data. However, for temperatures 
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Nf 





1 


2 


2+1 


3 


ro[MeV] 


270 


240 


208 


187 


178 



TABLE I: The critical temperature To for Nf massless fla- 
vors according to [s^. The value for 2 + 1 flavors has been 
estimated by using HTL/HDL theory for a massive strange 
quark with tUs — 150 MeV. 

T < 1.5Tq the deviations of the Polyakov loop potentials 
to the pressure are negligible. 

In the presence of dynamical quarks, the running cou- 
pling of QCD is modified by fermionic contributions. The 
size of this effect can be estimated within perturbation 
theory, see e.g. This leads to an 7V/-dcpendcnt mod- 
ification of the expansion coefficients in the polynomial 
Polyakov loop potential. The effect can be mapped onto 
an 7V/-dependent Tq resulting in Tab. H] [s^. Similar to 
[43| . the critical temperature Tq decreases for increasing 
Nf. 

In this work the To-dcpendence on the phase structure 
will also be investigated for the logarithmic and polyno- 
mial potential. 

III. THERMODYNAMIC POTENTIAL 

In order to investigate the phase structure of the 
PQM model the thermodynamic potential is evaluated 
in mean- field approximation similar to [36l . Isoj . In the 
derivation we will consider the dependence on the tem- 
perature and three quark chemical potentials, one for 
each flavor /. However, in the isospin-symmetric case the 
light masses are degenerated, i.e. m; = m„ = md and 



the thermodynamic potential depends only on two inde- 
pendent quark chemical potentials, = (/i„ + ^(i)/2 and 
fig. Note, that the quark chemical potential /x^ is simply 
related to the baryochemical potential fiB by = ^lb/^ 
with /is = 2^/ + /ij. Thus, only two order parameters 
(To and (Tg emerge in the singlet-octet basis. In the fol- 
lowing these condensates are rotated to the non-strange 
{ox) and strange (ctj,) basis, see [39] for details. 

In mean-field approximation the thermodynamic 
potential consists of the mcsonic U (ax, cry), the 
quark/antiquark flqq and the Polyakov loop contribu- 
tions 

n = U{ax,ay) + ngg{ax,ay,<^>,^) +U{<i>,^) . (18) 
Explicitly, the mesonic contribution reads 

Tll^ c 

U{ax,(Ty) ^ — {(tI + ay) - hxa-x - hyay - l^'^l^v 

+ ^^l^l + \ (2Ai + A2) a* + i (2Ai + 2X2) al , (19) 

and has six parameters w? ,c,Xi,\2,hx,hy. They 
are fitted to well-known pseudoscalar meson masses 
"^TT , m-K , rn^ + m^, and two weak-decay constants , Jk- 
The mass of the scalar a meson, rua, is used to complete 
the fit. Since the experimental situation concerning this 
resonance, cf. [i^l is not clear we vary 771^ from 500 MeV 
to 800 MeV which results in different parameter sets, see 
also [39} . 

The Polyakov loop variables are coupled to the 
fermionic part via 



n,,{ax,ay,<P,^) = -2T ^ y |ln [l + 3(a> + $ 



g-(B,,/-M/)/T^g-(E,,/--A'/)/T ^ g-3(B,,/-M/)/5 



f—u,d.,s 



In 



1 + 3('l + $e"(-^'-f+''-f)/'^)e"(^''-^+''-f)/'^ + e-^^^^-f+i'f''/'^' | . (20) 



The flavor-dependent single-particle energies are given by 

E,j = y^F + „,2 (21) 

with the flavor-dependent quark masses 

mi = gax/2 and = gayjspi (22) 

for the light and strange quarks, respectively. The 
Yukawa coupling g is fixed to reproduce a light con- 
stituent quark mass of « 300 MeV. This in turn 
predicts a strange constituent quark mass of TOs k, 433 
MeV. From Eq. (PH)) one recognizes the suppressions of 
single-quark contributions in the hadronic phase where 



$, $ ^ while baryon-like objects, which are composites 
of A^c = 3 quarks, contribute to the potential. In this way 
some confinement properties are mimicked in this model 
in a statistical sense. On the other hand, for $ ~ 1, the 
Polyakov loop modifications decouple from the fermionic 
part and Eq. (|20p turns to the usual quark/antiquark 
contribution of the pure quark-meson model [39| . 

Finally, the temperature- and quark chemical poten- 
tial dependence of the four order parameters for the chiral 
and deconfinement transition are determined as solutions 
of the corresponding equations of motion. These coupled 
equations are obtained by minimizing the grand poten- 
tial, Eq. (|18p . with respect to the four constant mean- 
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fields (a)^, {a)y, ($) and ($): 



an 

dcjx 



day 



on 



an 



, 



(23) 



where min = {a^ = {(Jx) , (Ty = {ay) ,$ = ($),$ = ($)) 
labels the global minimum. 



IV. FINITE-TEMPERATURE PHASE 
TRANSITIONS 

The (pscudo)critical temperatures for the various tran- 
sitions are determined as the inflection point of the cor- 
responding order parameters. The resulting critical tem- 
peratures at vanishing chemical potentials for the non- 
strange and strange chiral transitions as well as 
for the deconfinement transition Td for three different 
Polyakov loop potentials and various combinations of ma 
and To are summarized in Tab. |TT1 

In general, the inclusion of the Polyakov loop shifts the 
(non-strange) chiral transition at vanishing quark chem- 
ical potential, Hq = 0, to higher temperatures. In ad- 
dition, nia- and Tq affect the chiral and deconfinement 
transition temperatures and their interplay as follows: 
larger values of Too- and Tq shift again the transition to 
higher temperatures. A similar trend concerning the TOo-- 
sensitivity on the phase structure was observed in a pure 
three flavor quark- meson model [soj . 

For the chosen parameter sets we find > T^. This 
behavior is expected since the Polyakov loop suppresses 
fermionic contributions to the thermodynamic potential 
in the confined phase. However, these contributions drive 
the chiral phase transition and thus a delayed decon- 
finement transition always delays the chiral symmetry 
restoration. This becomes most obvious for small values 
of Too- 500 — 600 MeV, where the chiral transition oc- 
curs in the pure QM model already at < 150 MeV. 
As a consequence, the PQM model cannot reproduce the 
Wuppertal-Budapest scenario with T^<i^Td~T^. 

Only for the polynomial Polyakov loop potential with the 
extreme parameters Tq = 270 MeV and ra^ = 500 MeV 
the case Ty. < T^ is seen. However, the deconfinement 
transition is a very smooth crossover and overlaps with 
the chiral transition. 

For the logarithmic and Fukushima potential only few 
parameter sets (to^ = 500 MeV and 600 MeV) are found 
where both transitions coincide. Even in the deconfined 
phase the quark contribution to the thermodynamic po- 
tential is not sufficiently large to restore chiral symmetry 
immediately and thus > Td is obtained. 

In the following we choose a parameter set where both 
transitions coincide at /i^ = as suggested by the RBC- 
Bielefeld and HotQCD lattice calculations [3, Q- This 
is the case for = 600 MeV and Tq = 270 MeV 
where all three Polyakov loop potentials yield coincid- 
ing transitions. This choice also excludes modifica- 
tions caused by different values of TOo-. Moreover, re- 
cent studies [45| point to TO^-values which lie in the 
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184 


278 
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146 
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129 


238 






poly 


270 


800 


228 
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poly 


208 


800 
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297 


178 




poly 


187 


800 


197 


295 


161 




poly 


178 


800 


194 


273 


154 




poly 


270 


600 


204 


262 


204 
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poly 


208 


600 


179 


270 


179 




poly 


187 


600 


171 


267 


171 




poly 


178 


600 


166 


236 


166 




poly 


270 


500 


193 


254 


220 




poly 


208 


500 


168 


239 


168 


*- 


poly 


187 


500 


159 


237 


159 




poly 


178 


500 


156 


236 


156 




log 


270 


800 


228 


302 


208 




log 


208 


800 


208 


293 


166 




log 


187 


800 


204 


251 


150 




log 


270 


600 


207 


274 


207 
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log 


208 


600 


180 


286 


165 




log 


187 


600 


172 


283 


150 




log 


270 


500 


198 


264 


198 


* 


log 


zUo 


cc An 


ioo 


200 


ioo 




log 


187 


500 


159 


252 


150 




Fuku 




800 


214 


289 


193 




Fuku 




600 


189 


260 


189 




Fuku 




500 


179 


250 


179 


* 



TABLE II: The (pseudo)critical temperatures of the chiral 
transition for the light {T^), strange (T^) and the deconfine- 
ment transition {Td) a.t = for different rric and Polyakov 
loop potentials in units of MeV, evaluated for physical pion 
masses. In the last column the coincidence of the chiral and 
deconfinement transition is labeled with 



lower range of the broad experimental cr-resonance re- 
gion (400 MeV ... 1200 MeV) 0. 

For a proper comparison of the model results with the 
HotQCD d and RBC-Bielefeld lattice data BEl,!!!], we 
also adjust the pion masses accordingly. We denote the 
mass values of rriK = 503 MeV and to^t = 220 MeV which 
are used in lattice simulations as 'lattice point' and the 
measured values of rriK = 495 MeV and TO^r = 138 MeV 
as the 'physical point'. In the model calculation the tun- 
ing of the meson masses is achieved by adjusting the ex- 
plicit symmetry breaking parameters hx and hy accord- 
ingly. Heavier meson masses yield also slightly heavier 
constituent quark masses 

mi K. 322 MeV (300 MeV) , to, w 438 MeV(433MeV) 

where the values in parentheses are obtained at the phys- 
ical point. 

The increased meson masses shift the transition to 
higher temperatures, see Tab. IIIII The strongest shift is 
seen in the non-strange chiral transition since the change 
of the light quark masses is larger than that of the strange 
quark mass in this sector. Both, the non-strange and 
the deconfinement transition still coincide. The strange 
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FIG. 1: (color online) The Polyakov loop expectation value ($) (left panel) and the subtracted condensate Aj.s (right panel) as 
a function of temperature for different Polyakov loop potentials in comparison to lattice data of the HotQCD (p4 and asqtad, 
Nr = 8) [H and Wuppertal-Budapest collaborations (WB, Nr = 10, = 135 MeV) [I^. Solid lines correspond to larger pion 
and kaon masses as used in the HotQCD simulations, dashed lines to physical masses. 



quark sector is almost unaffected. 



In Fig. 1(a) the Polyakov loop expectation value (<f>) 
is shown as a function of temperature for three different 
Polyakov loop potentials. The solid lines correspond to 
heavier meson masses evaluated at the lattice point and 
the dashed lines to the physical ones. The expectation 
values are slightly larger for Fukushima's potential than 
for the logarithmic potential which exhibits a sharper 
crossover. The influence of the higher meson masses is 
mild. The observed shift may be explained by the larger 
T^/Tq ratio for heavier pion masses. The lattice data 
show a rather smooth transition with significant lower 
values for the Polyakov loop expectation value. For clar- 
ity we omit the error bars of the lattice data here and in 
the following figures. 

The sharper transition in the model as well as the 
higher values of the Polyakov loop are probably a rem- 
nant of the construction of the Polyakov loop poten- 
tials in the pure gauge sector, where a strong first order 
transition is obtained at T ^ 270 MeV. In the broken 
phase, the logarithmic potential comes closest to the lat- 
tice data, while the difference between the model results 
and the lattice is considerable in the restored phase for 
all potentials. 

Only for the polynomial potential the expectation 
value increases above unity which hampers its interpre- 
tation as an order parameter, related to the free energy 
of a static color charge. However, these findings are in 
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164 


250 






poly 


270 


600 


216 


263 


215 




log 


270 


600 


215 


275 


208 


* 


Fuku 




600 


200 


261 


196 


* 



TABLE 111: Pseudocritical temperatures similar to Tab. [TI] 
but for rrin = 220 MeV. See text for details. 



agreement with other model studies for three quark fla- 
vors, e.g. [U,!!^]. In a previous two-flavor PQM study 
($) > 1 was only observed at non-vanishing chemical po- 
tential [13 . For the other Polyakov loop potentials the 
logarithmic divergence avoids expectation values larger 
than one. 

In contrast to the Polyakov loop expectation value, a 
comparison of the chiral order parameter (qq) is more in- 
volved since some renormalization/normalization factors 
are unknown. Thus, a more suitable quantity, to compare 
with, is the ratio 



A, 



{ll){T)^im/ms){ss)iT) 
(ll){0)-{m/ms){ss){0) 



(24) 



which includes the non-strange (11) and strange (ss) con- 
densate. This corresponds to 



(T) - (K/hy) {ay) (T) 



(0) - ih^/hy) 



(0) 



(25) 



in our model calculation, where the quark condensates 
have been replaced by (ax) and {(jy) accordingly. Since 
the bare quark masses m; , rhg are not directly available, 
the ratio of the explicit symmetry breaking parameters 
hx^y is used instead. They are proportional to the bare 
quark masses up to a constant that drops out in the ratio, 
see e.g. [l|. 



In Fig. 1(b) the results for the subtracted condensate 
Ai^s are shown. The model calculations for Aj.s are in 
better agreement with the lattice data than the ones for 
the Polyakov loop expectation value. The lattice data ex- 
hibit a somewhat smoother transition than the model re- 
sults. The logarithmic Polyakov loop potential generates 
the sharpest transition while Fukushima's and the poly- 
nomial potential are closer to the lattice results. The in- 
fluence of the larger meson masses is more pronounced in 
the broken phase, where the transition becomes smoother 
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FIG. 2: (color online) The normalized pressure (left panel) and the interaction measure (right panel) as a function of tempera- 
ture. The model calculations (PQM model with various Polyakov loop potentials and the QM model) are compared to lattice 
data {Nt = 8, p4 and asqtad actions) from Q. Solid lines correspond to larger pion and kaon masses as used in the lattice 
simulations, dashed lines to physical masses. 



for larger meson masses. The data of the Wuppertal- 
Budapcst group suggest a much smoother transition. 
This effect is enhanced in the figure due to the transi- 
tion temperature normahzation. The transition temper- 
ature is much smaller in the Wuppertal-Budapest simu- 
lation than for the model or HotQCD simulations. Since 
the Wuppertal-Budapest data are obtained for physical 
quark masses they have to be compared to the dashed 
model curves which suggest an even sharper transition 
than the ones obtained for heavier pion masses. 



V. QCD THERMODYNAMICS 

In order to address the modifications of the Polyakov 
loop on the equilibrium thermodynamics we evaluate sev- 
eral thermodynamic quantities and compare the PQM 
model analysis with the QM model without the Polyakov 
loop and recent QCD lattice simulations. The bulk ther- 
modynamic observables such as pressure, energy and en- 
tropy density are sensitive to the change from hadronic 
to quark-gluon degrees of freedom and thus reflect some 
deconfining aspects of a QCD phase transition. 



Pressure 

The pressure is obtained directly from the grand po- 
tential via 



(26) 



with the vacuum normalization p(0, 0) ~ 0. In the model 
calculations it is the primary observable and other ther- 
modynamic quantities can be obtained from it by dif- 
ferentiation. Note that this procedure differs from the 
lattice approach, where for convenience the trace of the 
energy-momentum tensor is taken as the basis for the 



thermodynamic observables evaluation [8|. For example, 
the pressure is then obtained by integration over the trace 
of the energy-momentum tensor Q^^^ (46j 

p{T)^p{T,) + £ dT'^e^^ (27) 

with the integration constant p(Ti). In order to make 
the constant as small as possible the temperature Ti is 
adjusted accordingly, e.g., Ti ~ 100 MeV [431 p{Ti) 
is finally dropped. The suppression of this constant cor- 
responds to an approximation in the lattice data. 

Another normalization of the pressure lattice data con- 
sists in applying the pressure from the hadron resonance 
gas (HRG), PHRG(Ti), which provides a realistic descrip- 
tion of the pressure in the chirally brokenphase. 

Since pwaiTi = 100MeY)/Tf ~ 0.8 g this normal- 
ization shifts the lattice data to slightly higher values. 

In Fig. 2(a) the pressure for the PQM and QM models, 
normalized to the Stcfan-Boltzmann limit (SB) of the 
PQM model, is compared to lattice data. Explicitly, the 
pressure of — 1 massless gluonic and oi Nf massless 
fermionic degrees of freedom tends to 



NcNf , (28) 



at high temperature, which we use for normalization. Of 
course, no gluons are involved in the pure QM model 
and consequently, the first term in Eq. ([28| does not 
contribute. Therefore, a lower pressure is obtained in 
the QM model for high temperatures in the figure. The 
same conclusion holds for the PQM with Fukushima's 
potential. 

As expected, the QM model fails in describing the lat- 
tice data for all temperatures, while all PQM model re- 
sults are in reasonable agreement with the data. The best 
agreement around with the lattice data is achieved 
with Fukushima's potential but it fails for temperatures 
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FIG. 3: (color online) The energy- (left panel) and entropy density (right panel) similar to Fig. O 



above 1.5T^ as already discussed in Sec. Ill Al Note that 
for the HotQCD lattice data the transition temperature 
lies between « 185 - 195 MeV. We have used the 
average value = 190 MeV in the figures. Physical me- 
son masses reduce the pressure in particular around 



(dashed lines in Fig. 2(a) I. 

The better agreement of the pressure of Polyakov-loop 
augmented chiral models with lattice data has already 
been observed in two-flavor studies such as in [sB] for the 
PQM and in [l^l for the PNJL model. However, these 
model findings, obtained for physical pion masses, are 
compared with lattice data for two rather heavy quark 
fiavors p^ . To some extent, the agreement is there- 
fore fortuitous since, on the one hand, mesonic fluctu- 
ations are more suppressed in lattice simulations if heav- 
ier masses are applied and, on the other hand, mean-fleld 
approximations do not consider fluctuations around the 
mean-fields. In this work larger pion masses are taken 
into account for the lattice comparison. 



Interaction measure, energy density and entropy density 

The interaction measure A = e — 3p is related to the 
jrand potential via a temperature derivative 



dTKT'^ 



dT 



(29) 



and can be obtained directly in lattice simulations from 
the trace of the energy-momentum tensor 



9/^ 

rp4 



3p 



rp4 



A 

rp4 



(30) 



In the interaction measure finite volume and discretiza- 
tion effects are more apparent than, e.g., in the pressure 
(c.f. i). 

In Fig. |2(b)] the model calculations for the interaction 
measure are compared to ~ 8 lattice data. The height 
of the peak in the lattice simulations depends signifi- 
cantly on the lattice size as well as on the lattice action. 



Away from the phase transition the volume dependence 
is less pronounced [1]. In the chirally broken phase, the 
lattice data are close to all model curves. Around the 
transition at T « Fukushima's potential is closest to 
the asqtad-action data, while the remaining two model 
curves are in better agreement with the p4-action data. 
In general, the peak height decreases on larger lattices 
and in particular for the p4-action. For these tempera- 
tures the asqtad-action shows a weak iVT-dependence [1] . 
In contrast to the pressure, the logarithmic and polyno- 
mial Polyakov model version cannot describe the lattice 
data in the symmetric phase while Fukushima's ansatz 
approaches the data at least up to T ~ 1.5T^. 

In a PNJL study similar findings are obtained in the 
broken phase, but a more rapid decrease is seen in the 
restored phase [25j. For physical pion masses, the inter- 
action measure and the height of the maximum decreases 
for temperatures below the transition. Above the tran- 
sition the interaction measure increases for smaller pion 
masses. 

To complete the thermodynamic comparison we also 
show the energy density e and the entropy density s, de- 
fined as 



e = —p + Ts 



on 

df 



(31) 



in Fig.[3]as a function of temperature. The rapid increase 
of the energy density signals the liberation of light quark 
degrees of freedom. It increases almost to the value of an 
ideal gas of massless quarks and gluons. 



Speed of sound Cs and specific heat cy 

For hydrodynamical investigations of relativistic 
heavy-ion collisions the speed of sound Cs is an impor- 
tant quantity. Its square at constant entropy is defined 
by 



dp 




dp 
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(a) Specific iicat capacity cy jT^ 

FIG. 4: (color online) The specific heat capacity in the PQM 
model as a function of temperature. Solid lines have been 
obtained for m-n — 220 MeV, dashed lines for physical values. 
The arrows indicate the SB-limit. See text for details. 



where 



cv 



dT 



-T 



QT2 



(33) 



denotes the specific heat capacity at constant volume. At 
a second-order phase transition, where the specific heat 
diverges with the critical exponent a, the sound velocity 
drops to zero since the entropy density stays finite, see 
e.g. [i^ for a two-flavor QM model calculation. 

The specific heat calculated in the three-fiavor PQM 
model at = as a function of the temperature is 
displayed in Fig. 4(a) for three different Polyakov loop 
potentials. The solid lines in this figure correspond to 
masses used in lattice simulations and the dashed lines to 
physical masses. The specific heat grows with increasing 
temperatures, peaks at the transition temperature and 
approaches the corresponding SB-limit for high temper- 
atures. The SB limits deviate due to the differences in the 
contributing gluon degrees of freedom. For Fukushima's 
potentials the first term in Eq. ((28)) is neglected. For 
smaller masses the peak is sharper. Just after the peak a 
broad bump around 1.2 T /T^ is found for the polynomial 
and Fukushima's potential (for the latter with physical 
masses). Then cy goes gradually to the ideal gas value 
from above. 

The velocity of sound as well as the ratio p/e are shown 
in Fig.[5]as a function of temperature (left panel) and e^/"* 
(right panel). The latter plot eliminates uncertainties due 
to the normalization of the transition temperature. For 
all models including the pure QM model approaches 
the ideal gas value of 1/3 at high temperatures of about 
2.5Tc and can well be approximated by the ratio p/e. 
Furthermore, the value of p/e matches with the one of 
Cg for low temperatures. In between these two limits 
is always larger than p/ 1 except near T^, similar as in a 
two-flavor PNJL analysis [T5| . 

Near it drops to a minimum of about 0.04 for the 
PQM model, whereas the minimum is close to 0.16 if the 



Polyakov loop is neglected. While the latter value 0.16 is 
close to the results for pure YM theory on the lattice [4§| 
we find a minimum including the Polyakov loop almost 
half of the softest point P/e = 0.075 in [H^l for two-fiavor 
staggered fermions. It is interesting that in our work 
even at temperatures around O.Sr^ is not larger than 
0.1 if the Polyakov loop is included. This is in contrast 
to Ref. [5l| where a confinement model has been used 
and values of about ^0.2 around Q.bT^ and = 0.15 
around are found which correspond to our values for 
the pure QM model. 

Since the speed of sound is related to the interaction 
measure via 



A 

e 



ip 



I -id 



(34) 



a minimum in the speed of sound as expected near a 
phase transition translates to a maximum in the con- 
formal measure A/e. For high temperatures where 
reaches the SB limit of 1/3 the conformal measure should 
vanish. 

In the left panel of Fig. [S] where both quantities 
and p/e are shown as a function of e^/"* the best agree- 
ment with the lattice data is achieved with Fukushima's 
potential model. 



Quark number susceptibilities and Kurtosis 

Quark number susceptibilities incorporate information 
on thermal fluctuations of the degrees of freedom that 
carry a net number of light or strange quarks. They are 
also related to the second order cumulant of an expansion 
in the net quark number Sriq = Uq — {uq), e.g. j53 |. 

In the transition region these susceptibilities changes 
because the carriers of charge, strangeness or baryon 
number are heavy hadrons in the low temperature phase 
but much lighter quarks at high temperatures. Thus, 
at low temperatures they describe the fluctuations of 
hadrons carrying net light quark or strangeness quan- 
tum numbers. For example, the light quark number sus- 
ceptibility, xii receives basically only contributions from 
the lightest hadrons, the pions, at low temperatures, 
i.e. Xi/T"^ oc exp(— m^/T) where the index I refers to 
the light quark flavors. Correspondingly, Xs/T"^ receives 
contributions only from the lightest hadron that carries 
strangeness. For high temperatures and Nf quark fla- 
vors these susceptibilities approach the SB value for an 
ideal massless quark gas with Nc colors, i.e. Xq/T'^ ^ 

Furthermore, the quark number susceptibility Xq cor- 
responds to the first Taylor coefficient C2 of an expansion 
of the thermodynamic pressure in powers of [iq/T 



V_ 



OO 
Tl = 



.(T)(^)" withc„(r): 



1 d"{p/T*) 



n\ diiiq/TY 



(35) 
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FIG. 5: (color online) Left panel: The squared speed of sound (solid lines) and the ratio p/e (dashed lines) as a function of 
temperature. Right panel: For comparison with lattice data (p/e from 0|) both quantities are also shown as a function of e^''*. 
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(a) Light (solid) and strange (dotted) quark number 
susceptibilities xi,b 



(b) Ratio of light and strange quark number 
susceptibilities (2xs)/x! 



FIG. 6: (color online) Left panel: The light (solid lines/filled symbols) and strange (dotted lines/open symbols) quark number 
susceptibilities in comparison with lattice data (N^ ~ 8) with larger lattice pion and kaon masses @]. Fo r Y s we also show 
results of the Wuppertal-Budapest group (WB) obtained with physical pion masses on A^t ~ 10 lattices [ly]. Right panel: 
The ratio of the strange and light quark number susceptibilities with physical meson masses (dashed lines) and lattice masses 
(solid lines) in comparison with lattice data. Also shown are the results for a hadron resonance gas model (HRG) including 
resonances up to 1.5 GeV (upper curve) and 2.5 GeV (lower curve) from Q|. 



It can be obtained from the grand potential as follows 



2^ ^ 2c, - ^ 



(36) 



At a second order transition Xq diverges and, since the 
quark number density stays finite, the isothermal com- 
pressibility also diverges, see e.g. [53j . 

In the present (2-|-l)-flavor analysis we investigate the 
fight and strange quark number susceptibilities sepa- 
rately 



The required derivatives of the thermodynamic poten- 
tial can be calculated with a novel algorithmic differen- 
tiation technique, described in [s^ ]. 



Both susceptibilities, normalized according to their SB 
limit, are shown in Fig. |6(a) Around the chiral transi- 
tion temperature the light susceptibility xi raises faster 
than the strange susceptibility. In general, the PQM 
model results for the susceptibilities are below the lattice 
data. Closest to the model findings are the Wuppertal- 
Budapest data for Xs which have been obtained with 
physical pion and kaon masses. However, the influence 
of the pion mass in the strange sector is weak. The 
better agreement might be related to the larger lattices 
{Nr — 10) of the Wuppertal-Budapest group. 

To avoid unknown normalization factors in the lattice 
data (cf. [HI) it is more convenient to consider the ra- 
tio of the strange an d ligh t quark number susceptibili- 
ties as shown in Fig. |6(b)[ In the restored phase and 
around the transition the agreement with the lattice data 
is good. Note that the two different lattice actions devi- 
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FIG. 7: (color online) Left panel: Kurtosis i?, for uniform quark chemical potentials as a function of temperature. Right panel; 
Kurtosis Rs for the strange quark sector. The lattice data {Nt = 6) are taken from [3^ . 



ate strongly. 

The ratio of the strange and light suseeptibilities (right 
panel) shows strong deviations of the model calculation 
from either the lattice data or a HRG model calcula- 
tion in the broken phase. This might be attributed to 
the omission of mcsonic fluctuations in the mean-field 
approximation of the PQM model. Since the light and 
strange fluctuations arc mainly driven by the light pions 
and heavy kaons respectively, this results in rather small 
ratios for the HRG model calculation whereas values in 
between 0.5 — 1.0 arc found for the PQM model. 

Another quantity which probes deconfinement is the 
kurtosis of the quark number fluctuations Rq. It is de- 
flncd as the ratio of the fourth derivative of the grand 
potential with respect to ^q, i.e. Xq,4 = —d'^^/dfig, and 
the quark number susceptibility 



Xq 



(38) 



It is sensitive to the net quark content of baryon num- 
ber carrying effective degrees of freedom. At small tem- 
perature, the effective degrees of freedom in the PQM 
model are three-quark states since the Polyakov loop 
vanishes there and suppresses the one- and two-quark 
states in the thermodynamic potential, sec Eq. (|20| . 
Thus, the quark/ ant iquark contribution of the thermo- 
dynamic potential is basically that of a non-interacting 
system with particles and antiparticles carrying baryon 
number B = ±1. Consequently, the kurtosis tends in 
the low-temperature conflncd phase with \B\ = 1 to 
Rq = {^NcBY = 9. Without the Polyakov loop, e.g. in a 
pure QM model, the kurtosis would tend to one since in 
this case single quarks with \B\ — 1/3 would contribute 
at all temperatures. Thus, the correct low-temperature 
value for the kurtosis is obtained only with the Polyakov 
loop. 

At high temperatures the Polyakov loop tends to one. 
Fig. 1(a) As a consequence, the quark/antiquark con- 



the kurtosis reaches the SB value, Rq — )• G/tt^. Both 
temperature limits of Rq are precisely reproduced with 
the PQM model as shown in Fig. [T] The kurtosis is al- 
most temperature independent far below and above the 
transition. But around the transition it is strongly de- 
pendent on the meson masses and the effective imple- 
mentation of the Polyakov loop and its coupling to the 
quark sector of the model. The PQM model calcula- 
tion is consistent with lattice data {Nr = 6) [s^. At 
the transition the lattice data depict a peak behavior in 
Rq for larger meson masses (mjr ^ 220 MeV) but the 
volume dependence and errors there are still large. As 
already mentioned, only the PQM model provides the 
correct low-temperature limit Rq = 9 and reaches the 
high-temperature limit already at T ~ 1-3T^ similar to 
a two-flavor PQM analysis with a polynomial Polyakov 
loop potential j5a |. For Fukushima's potential we ob- 
serve a similar dependence of the peak structure on the 
pion mass, i.e., it vanishes for larger masses. The peak 
structure is more pronounced for a logarithmic potential 
and survives if the masses are increased. Interestingly, in 
a three-flavor PNJL analysis with Fukushim a.'s p otential 
no peak is seen even for physical pion masses [25| . A two- 
flavor PNJL model with a logarithmic potential gives a 
temperature dependent kurtosis in the broken phase for 
T < T^, which we also find in this work if we consider 
only the kurtosis for light flavors, i.e., Ri = T'^XiaIxi- 



tribution of the thermodynamic potential, Eq. (j20p . cor- 
responds to a system of non-interacting quark gas and 



In Fig. 7(b) the kurtosis for the pure strange sector 
Rs = T'^Xsa/Xs is presented. For the lattice simula- 
tions Nr = 6 and physical strange quark masses have 
been used. Note that chiral symmetry in the strange sec- 
tor is restored at ~ 1.3T^. For these temperatures 
the lattice data are already quite close to the SB-limit 
while the model analysis still yields an enhanced ratio. 
The strange quark kurtosis rises, and peaks at the chiral 
symmetry restoration for the light quarks. It stays above 
the lattice data until the symmetry is fully restored also 
in the strange sector and drops finally to its SB limit for 
high temperatures. 
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VI. SUMMARY 

In the present work we have extended a previously 
introduced two-flavor Polyakov-quark-meson (PQM) 
model [1^ to three quark flavors. This model with 
the Polyakov loop can address certain confinement is- 
sues and is an improvement of the known chiral quark- 
meson model with three quark flavors. Several aspects, 
in particular the bulk thermodynamics, of the (2 + 1)- 
flavor QCD phase transitions are in much better agree- 
ment with recent (2 -I- l)-flavor lattice simulations. 

The Polyakov loop potential is not uniquely deter- 
mined and different versions are available. Here, we have 
considered three different effective Polyakov loop poten- 
tials which all reproduce a flrst-order phase transition in 
the pure gauge sector of the theory at = 270 MeV. In 
two Polyakov potential versions a further parameter Tq 
enters which we have varied in order to investigate the 
influence of finite quark masses on the transition. The 
scalar and pseudoscalar chiral sector of the PQM model 
is fitted to the mass spectrum and decay constants in 
the vacuum. The input masses are experimentally well- 
known except for the scalar cr-resonance. Different sigma 
mass values have a strong influence on the chiral phase 
transition. In this work we have additionally considered 
various input values for the cr mass similar as in (39| . 

For some parameter combinations a coincidence of the 
chiral and deconflnement phase transition at vanishing 
chemical potential is found which we have used to study 
the bulk thermodynamics of the finite temperature tran- 
sition. 

The model results are confronted with recent lat- 
tice simulations for 2-1-1 quark fiavors of the RBC- 
Bielefeld [sl] and HotQCD Q collaborations. The latter 
lattice group uses lattices with temporal extent Nt = 8 
and the quark masses in these simulations are still above 
the expected physical ones. These quark masses corre- 
spond to a kaon mass of niK = 503 MeV and a pion 
mass of TO^ = 220 MeV. This difference in the masses 
influences the bulk thermodynamics and its effect is also 
addressed in this work. 

In addition, for the chiral order parameter and the 
strange quark number susceptibility lattice results for 
physical quark masses on Nr = 10 lattices of the 
Wuppertal-Budapest group [l^ are available and com- 
pared with the model analysis. 

We observe a nice matching with the lattice data for 
the chiral transition, although the transition is smoother 
on the lattice, in particularly for the Wuppertal- 
Budapest results. The lattice data for the Polyakov loop 
arc below the model data, especially in the deconfincd 
phase. 

A very good agreement of the lattice data with the 
model is found for the equation of state in the transition 
region and for temperatures up to l.ST^. 

Due to the omission of mesonic fluctuations in the 
mean-field approximation a deviation of the model and 



the lattice data is expected at low temperatures. Sim- 
ilarly, at high temperatures, the gluon dynamics is not 
fully addressed in the model which may also cause further 
deviation. 

The pressure is not very sensitive to different imple- 
mentations of the Polyakov loop potential and to details 
of the lattice simulations. This is in contrast to the in- 
teraction measure where Fukushima's potential describes 
best the lattice data in this temperature regime, espe- 
cially, with regard to the asqtad-action. 

For higher derivatives of the thermodynamic poten- 
tial such as the quark-number susceptibilities the picture 
changes. The model results are always below the lattice 
data. This is the case not only in the chirally broken 
phase which is understandable due to the omission of 
mesonic fiuctuations but also in the symmetric phase. 

Furthermore, there are also some normalization factors 
involved in the lattice data as well as some volume de- 
pendence. For example, the susceptibilities decrease from 
A^^ = 4 to A^^ 8 lattices. Interestingly, the Nr = 10 
data for Xs of the Wuppertal-Budapest group are very 
close to our model results. However, the normalization 
factors drop out in the ratio Xs/xi- ^ consequence the 
model results come closer to the lattice in the symmetric 
phase. In the broken phase the PQM model is far off the 
lattice and HRG results which might be an artifact of the 
mean-field approximation. 

In summary, we found that the PQM model gives a 
reasonable description of the bulk thermodynamics in the 
QCD transition region. The current implementation of 
the Polyakov loop into the chiral quark-meson model is 
certainly an improvement but still further developments 
arc required. It seems that the remnant of the flrst-order 
transition in the pure gauge sector modifies the chiral 
and deconfinemcnt transition at finite quark masses such 
that the transition becomes much sharper compared to 
lattice data. The sudden liberation of the fermionic de- 
grees of freedom in the model triggers this sharp behav- 
ior. This might change if one goes beyond the mean-field 
approximation. In this case the mesonic fluctuations also 
contribute to the thermodynamics and drive the phase 
transition. 
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